Predicting the statistics of wave transport through 
chaotic cavities by the Random Coupling Model: a 
review and recent progress 

Gabriele Gradoni, Jen-Hao Yeh, Bo Xiao, Thomas M. Antonsen, 
Steven M. Anlage, and Edward Ott 

Institute for Research in Electronics and Applied Physics, 
Department of Physics, 
Department of Electrical and Computer Engineering 
University of Maryland, 
College Park, MD-20742, USA 



Abstract 

In this review, a model (the Random Couphng Model) that gives a statisti- 
cal description of the coupling of radiation into and out of large enclosures 
through localized and/or distributed channels is presented. The Random 
Coupling Model combines both deterministic and statistical phenomena. The 
model makes use of wave chaos theory to extend the classical modal descrip- 
tion of the cavity fields in the presence of boundaries that lead to chaotic ray 
trajectories. The model is based on a clear separation between the univer- 
sal statistical behavior of the isolated chaotic system, and the deterministic 
coupling channel characteristics. Moreover, the ability of the random cou- 
pling model to describe interconnected cavities, aperture coupling, and the 
effects of short ray trajectories is discussed. A relation between the random 
coupling model and other formulations adopted in acoustics, optics, and sta- 
tistical electromagnetics, is examined. In particular, a rigorous analogy of 
the random coupling model with the Statistical Energy Analysis used in 
acoustics is presented. 

Keywords: Electromagnetic environment (EME), complex cavity, wave 
chaos, statistical modeling, impedance matrix, admittance matrix, 
short-range orbits, antenna, aperture, interconnected systems, statistical 
energy analysis (SEA). 
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1. Introduction 



The random coupling model (RCM) addresses the problem of statistically 
modeling the wave behavior of large irregular cavities connected to an exter- 
nal environment by one or more channels. The RCM is formulated in terms 
of impedance, a concept common to both the electromagnetic and acoustics 
communities. As such, many of the lessons learned in an electromagnetic 
context presented below have close analogs in the field of vibro-acoustics. 

Increasingly complex scenarios in electronics and telecommunications make 
the detail of structures and circuitry more and more difficult to model. In 
addition, often in optics, electronics, and acoustics, the need for higher bit- 
rates pushes wave sources to emit at very short wavelengths compared to 
the characteristics size of the excited system. In this regime, the scattering 
process can be very sensitive to details. A statistical approach then becomes 
appropriate. Specifically, one can ask what the statistics of quantities of in- 
terest are relative to a suitable random choice of the system. This is the goal 
of the RCM. Other statistical approaches have already been used to study 
reverberation chambers and wireless propagation channels. For example, a 
Gaussian state is assum ed a prior i for the complex Cartesian field upon 
maximization of entropy iHilll (120091 ) . 

The RCM approach involves two physical prescriptions inherited f rom the 
quant u m chaos theory ofcompound nuclei in nuclear reaction theory IWignei 
fll967h : lMitchell et all fl2010h . The first prescription is to replace the exact 
spectrum of the cavity with a spectrum whose statistics are obtained from 
those of a suitable ensemble of random matrices. The second prescription, 
motivated by the typical chaotic behavior of rays in complex enclosures, is 
to take for the mode amplitude a superposition of random plane-waves. 

The applicability of Random Matrix Theory (RMT) to the solution of 
simple wave equations in bounded regions w here the corresponding ray trajec- 
tories are chaotic was originally explored by iMcDonald and KaufmanI (119791 ) . 

Understanding the influence of ports (antennas/sensors) and apertures 
on cavity field dynamics is itself interesting in practice. The chaotic regime 
of fields inside an irregular cavity is established by the boundary geometry 
driving ray propagation. From experimental observations, it is known that 
this results in a very high variability of the wave properties of a system. The 
RCM incorporates this in a formally elegant and yet physically meaningful 
way. 

The main goal in formulating the RCM was to connect the behavior pre- 
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dieted by random matrix theory (RMT) developed for the case of isolated 
wave chaotic systems, e. g., closed cavities, to the study of the behavior of 
such a system when it is connected to other systems through deterministic 
wave propagation channels, and to investigate the resulting input/output sys- 
tem charact eristics. This wa s done in the quantum mechanical context, for 
example, by iBrouwerl (119951] . The RCM was developed for resonators with 



spatially localized channels 
for both the single channel 



Zhend fl2005 [) . In 2 006, studies were published 



Zheng et al.l (j2006al ). and the multiple channel 



Zheng et al.l (l2006b[ ) cases. Those two papers developed previou s studies co n- 
sidering isolated chaotic systems exhibiting Hamiltonian chaos lOtt (l2002h . 

The RCM attacks this problem from the point of view of port impedance/admittance 
matrices, in which the universal fluctuation of flelds within the cavity is joined 
to the specific radiation characteristics of deterministic channels, whose ge- 
ometry is presumed to be known a priori. 



2. Mathematical model 

The wave chaos inside the cavity is expressed as a universal (i.e., sys- 
tem independent) superposition of "chaotic" modes resulting in a universal 
ensemble of matrices, and the presence of channels connected with the exter- 
nal environment results in "dressing" this superposition with the radiation 
impedance matrix Z^"''^ of channel terminals. In the case of a cavity with a 
single port the cavity impedance Z'^"-'" can be expressed as 



(1) 



where Z^^""^^ is a port radiation impedance that will be defined later in the 
paper, is a statistically fiuctuating variable corresponding to a sum over 
chaotic modes, and Q {■} and 3fJ{-} denote imaginary and real parts. Thus, 
the statistics of Z'^"''" is determined by the statistics of the random matrix 
^ and by the nonstatistical matrix Z^"''^. The sum over chaotic modes is 
universal in the sense that it depends only on a single parameter c apturing 



Arnaut 



Wignei 



cavity losses and modal structure through the average quality factor 
( 2003 ) and the mean spacing between nearest neighbor eigenmodes 
( 119671 ) and on no other cavity details. 

Recently, the RCM has been developed for three dimensional cavities 
with distributed ports , e.g., complicated antennas, in the same form of ([1]) 



Antonsen et al.l (l201l[ ). The presence of apertures calls for a different mod- 



eling attack, leading to the same form of the RCM for the cavity admittance 
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matrix. The problem of external radiation coup ling; to the cavity through an 
aperture has been solved iGradoni et al.l (l2012bl Jal). The ROM introduces and 
incorporates system-specific characteristics efficiently. It will be also shown 
how deviations from pur e wave chaos can be included in the RCM by con- 
siderin g; short-ray orbits llshio and Burgdorferi (119951 ): iBaranger and Mello 
(119961 ). The RCM impedance expression ([T]) is preserved and the non- 
universal t e rm is augmented by tracing those "non-chaotic" short ray orbits 



Hart et all (12009f l. 



In realistic physical scenarios one may employ different kinds of channel 
coupling to a closed electromagnetic environment (EME). Imagine a situation 
involving compartments of ships/aircraft /automobiles/trains: there could be 
present very localized terminals belonging to circuitry located inside bays, as 
well as very large windows and ports. It is thus quite useful to formulate a 
"hybrid" RCM involving both port terminals and electrically wide apertures, 
and to consider interconnected cavities. 

In this paper, the general formulation of the RCM for ports and apertures 
is reviewed. Then, the predictions of probability distribution functions of de- 
sired input /output parameters (e.g., impedances, admittances or scattering 
parameters) obtained for single-mode channels by Monte Carlo computation 
of the RCM, and their experimental conffrmation by measurements on actual 
systems are recalled. Finally, an application of the model to a situation of 
practical interes t is discussed: the mod eling of linear chains of interconnected 
chaotic cavities iGradoni et al.l (j2012dl ). Interestingly, the formula obtained 
for the power ffowing through a chain of lossy chaotic cavities is similar to 
the statistical energy analysis (SEA) approach developed for the propagation 
of acoustic and vibrational waves. 



3. Historical remarks and methods 



Beginning in 1982 there was rapid development in the application of RMT 
to both quantum and classical chaotic systems. The ffrst use of RMT was 
devoted to reproducing the spacing distribution of energy levels in compound 
nuclei. Later Rydberg levels of the hydrogen atom in a strong magnetic ffeld, 
and elastomechanical eigenfrequencies of irregularly shaped quartz blocks 



were studied, as wel 
in a magnetic field 



as fiuctuations of c onduction through mesoscop i c wire s 



Guhr et al.l (Il998l ): iWeidenmiiller and Mitchell! (120091); 



Beenakker 


( 


1997 


); 


Alhassid ( 


2000) 



However, only 
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a very small number of experiments on wave-chaotic scattering and eigenfre- 
quencies of electromagnet ic cavities existed. The first s tudies of irregularl y 
shaped r aicrow a ye cavities IStockraann and SteinI (Il990l ) ; iDoron et al.l ( Il990l ) ; 
Sridhad f ll99lh : ICraf et all (Il992[ ) paved the way to electromagnetic wave- 



chaotic scattering research. Microwave cavities with irregular shapes (having 
ray trajectories evolved by a chaotic map) provided a simple and effective 
physical framework for the study of wave-chaos, where not only the magni- 
tude, but also the phase of scattering coefficients, can be directly measured 
from experiments. The picture of rays is analogous to that of particles bounc- 
ing in a confined system with irregular potential barriers, whose trajectories 
exhibit exponential divergence in phase space, i.e., two particles originating 
from the same position with slightly different linear momentum follow very 
different paths after a few bounces off the boundaries. This dynamic effect 
has an impact on the spectrum of complex wave systems that are chaotic 
in the classical limit, which can be modeled by an unpredictable (random) 
Hamiltonian. 

Th e derivation of the RCM explo i ts bo th the Wigner surmise Wignerl 
(119671 ) and the Berry hypothesis iBerryl (119771 ) . If the dimensions of the cavity 
are much greater than the excitation wavelength, the semiclassical regime can 
be invoked. Irregular boundaries create highly disordered mode topologies. 
Therefore, the mode amplitude is locally modeled as an isotropic random 
superposition of plane waves (the Berry hypothesis). The motivation of this 
hypothesis is that the orbits of rays in a fully chaotic system visit all regions 
of phase space ergodically. In particular, if some time along a long orbit is 
randomly chosen, the probability density of the orbit direction is uniform 
in to 271. At the same time, the resonant behavior of modes is preserved 
by using RMT to generate mode wavenumbers (the Wigner surmise). Here, 
eigenvalues of a random matrix of the Gaussian Orthogonal Ensemble (GOE) 
with Gaussian distributed elements have been used since the time-reversal 
invariance is assumed to hold. The present approach is not phenomenological 
as the invoked ensemble when using RMT e xpresses the sym metry properties 
of the chaotic class a cavity belongs to (see ISo et al.l (119951 ) for experimental 
work where "time reversibility" is violated by the presence of a magnetized 
ferrite). 

This motivates the use of the Gaussian ensembles. However, in the pres- 
ence of losses and intera ctions with the exter nal environment, e.g., reac- 
tions in compound nuclei iMitchell et al.l ( 120101 ). the universal properties of 



a system are intimately modified. For resonant electromagnetic systems. 
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a unified model describing the behavior of o p en systems with los s es wa s 
generated by several group s ISavin et al.l (120011): iFvodorov and SavinI (l2004f) : 

Fvodorov etaP fl2005h: ISavin et all fcoOSh. Sev- 



Savin and Sommers 



group s 

mm; 



Hemmadv et al 



eral ext er isions a n d experimental valid a tion were publishe d 

f|2005bl lal. l2006bl Jah: IZheng et al I fl2006ch : lYeh et al" (|2010bl lal.l2012b.a): lHemmadv et al 
(j2012|) related to the theorv of Izheng et all fl2006al lbllc]l After that, an effort 
was made to incorporate imperfections arising from the interacti on of ports 



with n earby objects and cavity walls, giving rise to short-orbits iHart et al. 
(120091 ) ■ More re cently, the mod e l has been extended to 3D cavities with 
distributed ports lAntonsen et al.l (120111). and to ac count for the interaction 
through arbitrary apertures iGradoni et al.l ( l2012cl ). Interestingly, the elec- 
trical formulation underlying the RCM is formally and physically similar 
to other formulations involving complex wave systems of different physical 
nature, e.g., acoustics. In the following, the analogy between the RCM and 
statistical approaches to understanding the properties of complex wave struc- 
tures in applied electromagnetics, acoustics & vibrations, and random lasers 
is briefly discussed. 

In statistical electromagnetics (SEM), statistics and stochastic analy- 
sis are used to model complex electromagnetic systems iHoUand and John 
( 1l999l ). The aim is, besides conceiving simple and effective design tools, to 
shed light on the collective behavior of modes and rays when subject to ir- 
regular and heterogeneous boundaries. Cavities and resonators constitute a 
natural and fertile terrain for SEM to be developed and applied. A num- 
ber of models have been proposed to explain and predict fluctuations inside 
cavities (see for example Hill ( 2009[ )). The RCM improves on those phe- 
nomenological the ories by describing random fluctuations through statistical 
ensemble theories Schwabl and Brewerl ( 2006 ). In contrast, the RCM builds 
on the exact model of impedance and admittance matrices as derived from 
Maxwell's equations. 

As previo usly remarked, the use of RMT in physical systems dates back 
to 1950, when lWigner and Dirad (Il950l ) adopted symmetry properties of sta- 
tistical ensembles to model the Hamiltonian of complex nuclear systems. Be- 
sides nuclear theo ry, early applications of RMT can be found in acoustics an d 
vibration systems ICouchraan et al.l ( 119921 ): iTanner and SondergaardI (120071 ): 
Wright and Weaverl ( 2010 ). The pure statistical perspective of physical sys- 
tems is even older and originates from thermodynamics. Interestingly, this 
perspective led spontaneously to the use of random matrices for modeling 
transmission of vibrational energy. Specifically, a model in acoustics that is 
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close to those in elect romagnetics can be found in the so-called statistical 
energy analysis (SEA) iLyonI ( 120031 ) . This method is currently used as a key 
strategy in developing nume rical tools for analysis of vibr ational energy in 
complex mechanical systems iMaksimov and Tannerl (120111 ) . 

Applications of RMT can be also found in the field of optics and pho- 
tonics. As an example, chaotic cavities are used to create several irregular 



Cao et al. 


(1999 


); 


Tiireci et al. 



(l2006h . 



4. Theory 



The RCM is conveniently formulated in terms of an impedance or admit- 
tance matrix. The use of network theory for s olving electro r nagne tic problems 
is widely accepted in engineering and physics iFelsen et al.l ( 120071 ) . This often 
leads to theoretical formulae that can be efficiently computed and compared 
with actual measurements. In particular, the use of open-circuit parameters 
such as impedances and admittances does not require any specific knowledge 
of the sources and detectors connected to the terminals where voltages and 
currents are supposed to exist. Imagine injecting an excitation current into 
a small antenna radiating inside an irregular cavity. The multiple refiections 
and scattering of rays off the metallic (lossy) walls establish a reaction to the 
current density flowing on the antenna surface. This creates an electromotive 
force that builds up the response voltage at the terminals of the antenna. In 
the case of a single port, the ratio of the reaction voltage over the excita- 
tion current deflnes the terminal impedance. This parameter has affinities 
with the concept of acoustic impedance, defined by the ratio of the acous- 
tic pressure, analogous to the reaction voltage in electromagnetic systems, 
and the acoustic volume flow, analogous to the electrical excitation current. 
The acoustic impedance at a particular frequency indicates how much sound 
pressure is generated by a given air vibration at that frequency. 

In the specific case where t he antenna ra d iated in side a cavity, the impedance 
is called the cavity impedance IZheng et al.l (j2006al ). In the presence of mul- 
tiple ports, the response effect of the excitation current on each antenna 
terminal can be ea sily separated by defining a cavity impedance matrix 
Zheng et al.l (l2006al ). 



(2) 
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To find in 3D cavities, a solution of the vector wave equation satisfying 
cavity boundary conditions is constructed. This is done through the usual 
procedure of expanding the cavity fields in a basis of electromagnetic and 
electrostatic eigenmodes. At this point, the prescriptions of the RCM based 
on wave chaos and random matrix theory are used. Firstly, the exact cavity 
spectrum is replaced with a spectrum of eigenvalues generated by a large 
random matrix of the GOE. The symmetry properties of the adopted matrix 
refiect the physical properties of the actual system, whence other ensembles 
can be invoked for generating the random eigenvalues, e.g., the Gaussian 
Unitary Ensemble (GUE) in case the time-reversal invariance is broken by 
the presence of a magnetized ferrite. Secondly, the mode amplitude (which 
is unknown) is replaced with a superposition of random plane-waves. Those 
two ingredients define what we call a "chaotic mode" of the irregular cavity. 

Starting from a deterministic mode expansion, the RCM prescriptions 
lead to a fiuctuating impedance matrix having dimension NpX Np, where Np 
is the number of ports connected to the cavity, viz., 

g^^^ = l^™"!} + [3? 1^™"'}] • ^ ■ [3? 1^""'^}] , (3) 

where Z'^"'^ = RJ''^'^ + (Z''"'^) is in the simplest theory an Np x Np diagonal 
matrix whose elements are the complex radiation impedances of each port. 
Here, the radiation impedance provides the linear relation between voltages 
and currents at a port in the case in which waves are allowed to enter the 
enclosure through the port but not return, as if they were absorbed in the 
enclosure. The matrix ^ naturally incorporates cavity homogeneous losses, 
and it is defined as ^ 

where /C^.-^ = k^,^/Ak'^, and a is the loss parameter defined below. In 
the lossless case, ^ is found to be an element of the Lorentzian ensemble 



Brouwer and Beenakkerl (1l997l ). Here, is a vector of uncorrelated, zero 
mean, unit width Gaussian rand om var i ables, and k"^ are the eigenvalues of 



a matrix selected from the GOE lMehta] (Il99ll ). where the central eigenvalue 



is shifted to be close to fcg = w^/c^, and u is the frequency of excitation. 
In case of high number of overlapping modes, e.g., occurring in overmoded 
reverberation chambers, the eigenvalues have unifor m distribution and are 



distributed as per Wigner's "semicircle law" (see ( iHemmady et al.l . 12012 
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Appendix) and iMehtal (jl99ll ) for more details). Eigenvalues associated with 
participating modes are scaled so that the average spacing between eigenval- 
ues near the central one is Afc^, which is selected to match the mean spacing 
of resonances of the enclosure in the frequency range of interest. The effect of 
wall losses, internal homogeneous losses and the cavity dimensions (volume), 
are embedded into a single dimensionless loss parameter 

appearing in which is formally identical to the inverse finesse parameter, 
widely used in optical resonators to quantify th e gran ularity of resonances 



over a frequency band of interest iBorn and Wolj (119991 ). and where the (av- 



erage) cavity quality factor is defined as 

Q = ^^. (6) 

with E the average energy stored inside, and Pd total power dissipated 
throughout the cavity. The greater a the higher the resonance overlapping 
over the considered bandwidth. 

It is through the matrix ^ that the propagation of waves in the enclo- 
sure from one port to another and back is modeled. Such a propagation 
is universally fluctuating with statistics regulated by a single parameter, a. 
The RCM treats ports by means of a slowly-varying-in-frequency radiation 
impedance matrix. The enclosure itself is modeled by (jl]) that has the phys- 
ical meaning of a normalized impedance, modeling the chaotic scattering 
taking place throughout the cavity. Rapid variations of impedance with 
frequency arise from this chaotic behavior. Similarly, in acoustics and vi- 
brational systems, the acoustic impedance usually varies strongly when the 
frequency is changed. The form ([3]) makes a perfect link with well-established 
deterministic theories. The novelty of the RCM resides in the definition of a 
'random' impedance that is dressed with deterministic features. 

In practical mode-stirred chambers, it is expected that perfect mixing 
leads to a zero mean complex field. However, an offset in the scatter plot of 
the measured field is often experienced. By adopting the Poisson kernel the- 
ory, the RCM has been extended to account for short-range ray trajectories 
(short-orbits) between the ports. Those short orbits create system-specific 
properties resulting in nonuniversal frequency fluctuations, and can be ac- 
conted for with an extension of the RMT. 
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The philosophy of the Poisson kernel Mello et all (llQSSi l ; iKuhl et all (120051 ) 
is to include those scattering matrix contributions arising from parts of the 
system that are non-random. This same perspec tive has be e n app lied to the 
impedance matrix Z of a microwave enclosure iHart et al.l (|2009[ ). Experi- 
ments involving mode-stirred cavities and ray-chaotic cavities often rely on 
the variation of boundary configuration to create several s tatistically indepen- 
dent r e alizations of the syst e m and to compile s tatistics iKostas and Boverie 
(I1991I ): ISchafer et al.l (120051 ): iKober et al.l (120111 ). However, i mperfections in 



the re alization of perfectly diffused fields have been observed iPrimiani et al. 



(120091 ) numerically and experimentally. This resulted in the beginning of 



a serious effort to m o del a n d understanding t hese "unstirred" occurrences 



Primiani and Moglid (120101 ): iPirkl et al.l (l201lh . From a wave chaotic per- 



spective, it is noted that the channel- channel interaction introduces deter- 
ministic field components that remain fixed throughout the ensemble (such as 
wall proximity on scattering objec ts inside the enclosure) desp ite the design 



of a very efficient stirring structure iMoglie and Primianil (|2012|). This kind o f 



"imperfection" has be en widely investigated in quantum chaos iPrangd (120051 ); 



Bulgakov et al.l ( 120061 ). Therefore, relevant ray trajectories, which remain un- 



another system-specific feature. Hart et al. 



Hart et al. ( 


2009 


); 


Yeh et al. 



leave a port and soon return to it, or another port, instead of ergodically 
sampling the system. These ray trajectories are also called "short orbits" 
or "short trajectories" . It is imagined that each ray leaves from a port p at 
Ti with a certain kj vector and it arrives at a final position r/ with a wave 
vector k/. It is further assumed that the ports are separated from the wall 
much more than a wavelength. 

The impedance form of the RCM is slightly modified by the presence of 
short-ra nge trajecto r ies. I n particular, it is found by a semiclassical theory 



of orbits iHart et al.l (l2009f l 



(7) 



where T is the Bogolomny transf er operator which in the electromagnetic 
case is given by iBogomolnyl (1l992l ) 



^fe,g/,k) = — 



n piS{ri,ri,]i)-ii 



CO?, Of 
cos Oj 



(8) 
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where 0i{6f) is the angle between the initial (final) wave vector and the 
surface at the position it leaves (hits), S" (r,, r^, k) is the classic action along 
the direct trajectory from to r/, and ^Dj.-j.j is the stability of the orbit 
from Vi to Yf (monodromy) that physically represents a geometrical factor of 
the trajectory taking account of the spreading of the ray tube along its path. 
It is also found that 

'g^ ■U = ^{g""^] . (9) 

Ultimately, the correction introduced by the presence of short-orbits in the 
RCM average impedance as 



.mdl 1/2 



where the elements of z are expressed as 



■ad^ 1/2 



(10) 



b{n,m) 



Ikl + A; L 



IklL 



port 



:ii) 



with b (n, m) an index over all classical trajectories which leave the n port, 
bounce (3b times, and return to the m*^ port. Lb is the length of the trajectory 
b, k = k/ {2Q) is the effective attenuation parameter taking account of 
loss through the average quality factor Q, and Lport is the port- dependent 
constant length between the n*^ and m*'* port. The geometrical factor Db is 
now a function of the length of each segment of the trajectory, the angle of 
incidence of each bounce, and the radius of curvature of each wall encountered 
in that trajectory. The application of this theory to reverberation chambers 
could pave the way to understand and use the single-frequency regime, where 
comparison of experimental and theore tical statistics often e xhibits strong 
deviation from asymptotic distributions iPrimiani et al.l (|2009[ ). 

Researchers have examined short orbits in cases where t he syste r n and 
the ports can be treated in the semiclassical approximation iPrangd (120051 ) 
or considered the effect on eigenfunction correla t ions d ue to short orbits as- 
sociated with nearby walls lUrbina and Richterl (l2006l ). Theoretical results 
have been previously obtained th at acknowledge the effects of short-orbits 
in the context of quantum graphs iKottos and SmilanskyI (I20031). Th e RCM 
for micro^yave cav ities has been confirmed by experiments iHart et al.l (120091 ); 
Yeh et al.l (l2010al ). Moreover, experiments in a 2D microwave cavity have 
extracted a measure of the microwave power that is emitted at a certain 
point in the cavity and returns to the same point after following all possible 
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Figure 1: Wave chaotic enclosure with both one aperture and one port. 



classical trajectories of a given length IStein and StockmannI ( 119921 ) . None 
of this prior work developed a general first-principles deterministic approach 
to experimentally analyze the effect of short ray trajectories. The effect of 
short traj e ctories in two-port wave-chaotic cavities has been demonstrated 
Yeh et al. (l2010bl Jal) . The results of the short-orbit correction can be general- 
ized to distributed ports and apertures by similar procedures. For apertures, 
it is interesting to think about the proximity of a dielectric/metallic object 
that is immersed in the wave-chaotic field. This physical situation has been 
investigated by Harrington in the case of regular cavities , where the a pertu re- 
object resonance was predicted by method of moments iHarringtonl (119821 ). 



4-1- Aperture radiation of cavities 

Besides localized ports, the RCM can also include a description of aper- 
tures, a situation depicted in Fig. [H Ongoing research is focusing on the 
extension to very large and irregular apertures, as well as to very small aper- 
tures, where singular fields strongly affect their (system specific) radiation 
Bethel (jl944j ). The aperture excitation of an irregular enclosure is indeed a 
rather new theoretical problem to investigate. Interestingly, such a theory de- 
velops well on the existing literature on cavity backed apertures, and it links 
deterministic theories with statistical theories on reverberation chambers. 
The work of Harrington on cavities showed that boundary-value problems 
involvii ig apertures are conveniently described in terms of the admittance 



matrix 



Harrington and Mautzl (jl976[ ). Typically, radiation admittances are 



obtained by numerical computation of cavity eigenmodes. On the other hand. 
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it is known that when the cavity is strongly perturbed by irregular objects 
and/or boundaries, modes change the ir topology and spectru m completely, 
mode stirred enclosures/chambers iBunting and Yu| (12004 ). and the fields 



e.g 

become amenable to statistical description Hill (2009i) . Also in the presence 
of apertures, thanks to the RCM perspective of separating system-specific 
from universal characteristics, those two apparently distinct approaches can 
be unified through wave chaos theory. Similar to a port, the aperture geom- 
etry is fully included in the model through the radiation admittance matrix. 
The procedure to be followed is practically the same that has been used 



to arrive at Eq 



(13]). A detailed mathematical derivation is described in 



Gradoni et al.l (j2012d ). and is briefiy recapitulated below. 

The aperture has been treated as an opening in a zero thickness metallic 
plane, and the analysis started by expanding the aperture (tangential) field 
distribution in a basis of modes 

Et = J2Vses{^±) , (12) 

s 

where Vs are aperture voltages representing the electric mode amplitudes, 
is the spatial coordinate within the aperture the aperture modes 

and have only transverse fields, normalized such that f , dx"^, leoP = 1. 

' J aperture -L i 

Solving Maxwell's equations with Dirichlet boundary conditions = 0, 
yields an expression for transverse magnetic fields in terms of the same basis 
of modes 

H( = ^ Isfi X e, (xx) , (13) 

s 

where fi is the outward normal to the cavity, which has been taken to be in the 
2;-direction. The linear relation between magnetic field mode amplitude Ig, 
and electric field mode amplitude Vg, is expressed in terms of an admittance 



(14) 



Here Y^'} = y'-™'^^ or Y^'} = Y^T^^ depending on whether the aperture ra- 

ss ss ss ss 

diates into free space or into a cavity. In particular, in presence of the 
cavity, the boundary-value problem is solved by expanding the cavity fields 
in modes, and then by applying the RCM prescriptions for the properties of 
these modes. In this case, the RCM formulation leads to a specific relation 
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Figure 2: Geometry of an aperture of arbitrary shaper illuminated by an 
external plane-wave. 



between free-space admittance Y^^i"''^^ and cavity admittance Y^T'"\ viz., 



(15) 



where the matrix ^ is the same as defined in (jl]), and the elements of Y''"-'^ 



are defined by iGradoni et al.l (l2012d ) 



4-2. Coupling with external radiation 

In practical electromagnetic systems, an aperture can be fed by a waveg- 
uide with the same transverse sectional shape, or it can be exposed to ex- 
ternal radiation. In both cases the tangential field expansions ( !T2|) and (fT3!) 
can be preserved. It is thus instructive to analyze the typical scenario de- 
picted in Fig. [2] of an oblique plane-wave with wave vector k*"*^ and polar- 
iz ation of magnetic field h'"''^ (p erpendicular to k*"'^) exciting the aperture. 



In iHarrington and Mauta (Il976l ). this analysis has been performed by con- 
tinuity of tangential magnetic field at the aperture plane, that leads to the 
equivalent network equation by projecting the two magnetic field express ions 
on the aperture basis, and equating the amplitudes Gradoni et al. ( 2012b al, 3) 



(16) 



where V is a vector containing the voltages Vg, PJ^^ = —n ■ (— k^'^) x h*"^, 
and e is the Fourier transform of the aperture basis mode for the electric 
field. 
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4-3. Hybrid formulation: ports and apertures 

The use and verification of tlie RCM for apertures tlirougli quadratic 
(power) measurements requires tlie presence of antennas/ports inside tfie 
cavity. Tlierefore, it becomes natural to formulate a hybrid RCM including 
both the aperture admittance matrix and the port impedance matrix. 

The RCM has been extended to account for the joint presence of elec- 
trically wide openings of arbitrary shape and localized ports. A schematic 
framework of such a situation is reported in Fig. [H In that case, an input 
column vector has been constructed that consists of the aperture voltages 
and terminal currents, and an output vector ip consisting of the aperture 
currents and terminal voltages 



and 



Lp 

La 
Vp 



(17) 



;i8) 



where V^,p are the aperture and port voltages, and La,p are the aperture 
and port currents. These are then related by a hybrid matrix LL, = LL ' 4*^ 
where 

1/2 / X -1 1/2 



Here the matrices U and V are block diagonal, viz.. 



(19) 



U 



Yrad g 



(20) 



and 



Z = 3^(£) ■ (21) 

The dimension of f/ and V_ is (A^^ + Np) x (A^^ + Ap), where Np is the num- 
ber of port currents and A^, is the number of aperture voltages. The linear 
problem at hand is thus formulated in terms of matrices with relatively small 
dimensions, typically given by the number of ports and aperture modes used 
to access the electromagnetic structure. Stri ctly speaking , the RCM formu- 



lation has facilitated a model order reductio mYelseji et al.l (120071 ) if compared 
to full-wave or other statistical models of cavity fields. Usually the superpo- 
sition of plane-waves or modes requires the computation of the field at each 
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point of the structure, thus producing very large matrices. In the different 
perspective of looking at channels, the available number of degrees of free- 
dom (typically constituted by the number of excited modes of the cavity) is 
reduced to the number of port profiles/ aperture basis functions. 



5. Random coupling model in action: numerical simulations 



Ensembles of cavities ha ve been gener a ted by the application of the Monte 
Carlo method described in IZheng et al. f l2006al Jbl): iHemmadv eFaP f l2012h . 
Once the system-specific properties of channels are predicted or measured 
separately under an a — )■ oo (free-space) condition, Monte Carlo simulations 
of (jl]) allow for generating an ensemble of cavity impedances of the form 
03]) or, identically, of cavity admittances of the form (ITS]) . A proper num- 
ber of realizations must be chosen to create adequate statistics close to the 
asymptotic behavior, which is Lorentzian in the lossless regime (a = 0), and 
Gaussian in the high-loss regime (a ^ 1), see Fig. (j3]). First of all, the 
universally fluctuating impedances (j4]) is recast as 



1 



T 



TT — A — ial — 



(22) 



The matrix $ is a M x iV coupling matrix with each element representing 
the coupling between the i^^ port profile/aperture mode {1 < i < M) and 
the j^^ chaotic mode of the cavity (1 < j < A^). Each $jj is an independent 
Gaussian-distributed random variable of zero mean and unit variance. The 
matrix is the transpose of and 1 is a x identity matrix. The 
matrix A is an A^ x A^ diagonal matrix with a set of N random eig envalues 
based on the GOE nearest-neighbor (Wigner) spacing distribution IWigner 
(119671 ). By repeating this procedure many times in order to give a suffi- 
ciently large ensemble of Z, the statistical description of the cavity is thus 
generated. Conversely, starting from cavity measurements of the impedance, 
the universal fluctuation of the single port impedance can be calculated by 
normalization of ex perimental result s , to com pare to RMT predictions of the 
universal behavior Hemmady et al. ( 2005bl a). Fig. |3] reports the distribu- 
tion of the real (a) and the imaginary (b) part of the normalized impedance 
parametrized as a function of the loss factor a. The loss parameter of the 
cavity at hand is predicted from the Weyl law on the basis of the volume 
of the cavity, its frequency of excitation {u = kc), and the electromagnetic 
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Figure 3: Random matrix theory prediction for the real (a) and the imaginary 



(b) part of the normahzed impedance z 
a. 



for various values of 



losses. For a 3D enclosure it is found 



a 



(23) 



where V is the cavity volume, and Q the average resonator quality factor. 

The effect of increasing the loss parameter results in reducing the magni- 
tude of fluctuation of impedance elements. It is to be noticed that for a = oo 
we retrieve the case where ports/apertures radiate in free-space. In the limit 
of infinitesimal losses the impedance fluctuation become very strong: t his 



is a situation that can occur in superco nducting cav i ties lAlt et al.l ( 119981 ). 



Tiireci et al.l fcoOGh 



m microwave 



in low- loss optical systems, e.g., lasers 
cavities ope rated in the non-Ericsson (weak overlapping) regime iDietz et al 
( 120081 . 12OO9I ) , and in reverberati o n chambers /e r iclosu r es operated in t he un - 
dermoded regime Arnaut ( 2001 ): Warne et al.l (2003); Orjubin et al. ( 2006 ): 
Arnaut and Gradonil ( 120111 ). 

In the limit of high losses, the diagonal (off-diagonal) elements of ^ become 

unit (zero) mean Gaussian random variables. 

Figs, m (red curves) reports a few experimental results for the frequency 
dependence of the real (a) and imaginary (b) single-port impedance measured 
at the coaxial port feeding a "bowtie" ray-chaotic cavity. A high frequency 
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Frequency (GHz) Frequency (GHz) 

(a) (b) 

Figure 4: Resistance (a) and reactance (b) of the measured port (input) 
impedance in a "bowtie" cavity. The sohd red curves are measured data. 
The dashed blue curves are created by combining the short orbit effect and 
the radiation impedance (thick black curves). 



variability due to fluctuations that ride on top of several systematic trends 
can be noticed. Rapid fluctuations come from longer orbits. The first sys- 
tematic trend is shown as a black solid line and represents the measured 
radiation impedance of the port. The second systematic trend is shown as 
the thin dashed curves which are semi-numerical results that contain the 
radiation impedance (measured) and the effect of short (less than 200 cm) 
orbits. This should be compared to the characteristic dimension of this bil- 
liard, which is £ = ^J~A ^ 34 cm. Note that the short-orbit curve follows the 
major trend in the data, but does not describe the rapidly varying part of 
the impedance. By removing all systematic effects up to an equivalent short 
orbit length of Cm = 200 cm one can obtain a good corr espondence betwe en 
the theoretical and experimental normalized impedance lYeh et al.l ( l2010al ). 

The question now arises: how far should the short-orbit calculation be 
carried? At short times only a few of the shortest orbits will contribute, 
whereas the number of orbits contributing to the impedance will grow expo- 
nentially as the time increases. The relevant time scale for terminating the 
short-orbit calculation is the Ehrenfest time. The Ehre nfest time is a time 
scale associated with the classical-to-quantum crossover lAleiner and Larkin 
(119961 ). and it is the time when wave packets start to deviate from the mo- 
tion of deterministic classical ray trajectories to fully-developed wave chaos 
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Figure 5: Network model of a chain of coupled cavities: the first cavity is 
excited by a localized source (port 1), and the coupled field at the last cavity 
is picked up by a localized detector (port N). The cavity-to-cavity coupling 
can take place through either single-mode channels (transmission lines) or 
multimode channels (distributed ports or apertures). 
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Schomerus and JacquodI (l2005l ). Therefore, it is meaningless to compute the 
orbits whose lengths correspond to time scales much longer than the Ehren- 
fest time. For a microwave cavity, the Ehrenfest time can be calculated as 

<.4.n(f), (24) 

where h is the largest Lyapunov exponent of the classical ray dynamics in the 
cavity, C is the characteristic dimension of the cavity, and A is the wavelength 
of the probing waves. The Ehrenfest times of the bowtie cavity and the cut- 
circle cavity have been estimated. The Lyapunov exponent for the bowtie 
billiard is h = 0.7737 ns~^, so the Ehrenfest time at / = 10 GHz is 3.2 
ns which corresponds to a length of 94 cm in free space. For the cut-circle 
billiard we found h = 0.6411 ns~^, so the Ehrenfest time at / = 10 GHz is 3 
ns which corresponds to a length of 92 cm in free space. 



6. Interconnection of cavities 



The RCM has been pushed further ahead to address the interconnection 
of complicated sy stems. Hence, t he scen ario of a linear chain of cavities has 
been investigated iGradoni et al.l (l2012dl ). 

The analyzed physical framework consists of interconnected two-port 
cavities, as shown in Fig. [51 A fundamental quantity of interest for this 
analysis is the input impedance of the chain. This can be found in an iterated 
manner starting at the last cavity, where Z^^^ = z[i \ to the first cavity, to 
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find 

The concept of trans-impedance has been adopted as a basehne to analyze 
and quantify the amount of signal coupled between two or more arbitrary 
elements (subsystems) of the chain. In particular, for the n-th cavity, the 
voltage at the first port of the n + 1-st cavity V^"^^^^ is related to the current 
exciting the first port of the n-th cavity as 



4"^ = \t- = 7T ■ (26) 

J j(n) y{n) ^(n+l) ^ 



,{n+l) y{n+l) ry(n) 

1 ^ ^in ^; 

+ zi 

Consequently, the ratio of power entering the A^-th cavity to that entering 
the first cavity can be expressed product, viz., 

p{N) 

= ~77vr ' (2'^) 



P. 



(1) 



where each impedance involved in ( l26l) is expressed by the RCM. 

The case of statistically identical cavities having moderate/high losses 
= «2 = • • • = «Af = > 1, deserves special attention. Correspondingly, 
the weak fiuctuation approximation can be employed leading to a simplified 
expression of (126|) . More precisely, in the high- loss regime, zf^-* ^ , and 
fiuctuating impedances can be tackled with first-order perturbation theory. 
This procedure results in a simplification for the denominator of ( l26ll . and 
leads to an expression for the coupled power ratio Rn that factorizes 

f{n,n+l) 

(") - n -^"^ ' (28) 

n=l 

where x*-"^ are exponentially distributed random variabes 

fx (x) = exp i-x) . (29) 
The transmission factors in (1281) are given by 



p(") p("-+i) 

' 7(") 1 7('T'+1) ^ 
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0.9 p > > > > 

8 - " = 6 N=2 „^ 




Figure 6: Probability density function of power ratios -Rat (a = 6) = 

for chains of up to seven cavities: the high-loss regime is assumed for all the 
(statistically identical) cavities in the chain. 



modeling the coupling between the n-th and n + 1-st cavity, and the loss fac- 
tor is given by ( 123|) . In iGradoni et al.l (l2012dl ). the Monte Carlo method de- 
scribed in section |5] has been used to generate ensembles of impedance matri- 
ces and compute the coupling in chains of cavities. Basically, the impedance 
matrix for each single cavity has been generated, and used in the theoreti- 
cal expressions ([3]) and fl22l) . Then, the power ratio has been evaluated for 
sequences of cavity chains of different length. For these studies, the cavities 
were assumed to be statistically identical. The PDFs of the logarithm of the 
power ratio for the case a = 6 is shown in Fig. Ei The high loss case (Fig. [H]) 
is well approximated by the analytic formulas fl25]) . fl2^ and fl5U]) . The distri- 
bution approaches log-normal as the number of cavities becomes large. This 
is evident from the comparison of the 7-cavity chain with a best-fit log-normal 
distribution. This behavior is expected since according to (|28l) and ( !29l) the 
power ratio becomes a product of independent and identically distributed 
random variables, and thus, as N becomes large the logarithm of the power 
ratio will be normally distributed. Anderson localization in the lossless case 
has been observed with respect to the number of ch ain elements by adoptin g 
a cascaded RCM in a chaotic dynamics perspective iGradoni et al.l (l2012dl ). 
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6.1. Connection with acoustics and vibrations 

It turns out that the flow of electromagnetic energy can be expressed in 
a similar way to that used to model the flow of vibrational wave energy. The 
global system is divided into complex subsystems intercon nected with each 



other through simple deterministic structures. In particular, iLyonI (|2003[ ) ap- 
plied a "thermodynamical" approach to arrive at a statistical model of energy 
exchange between two parts (subsystems) of the system, namely statistical 
energy analysis (SEA). Specifically, by using energy equipartition arguments, 
the average power passing from subsystem i to subsystem j (Pij) can be 



expressed as (ITanner and Sondergaardl . 120071 . Eq. (92)) 



P., = ud^,, (^-^) > (31) 
V di dj ) 

where uj is the mean frequency of the source, r]ij is the coupling loss factor, 
is the mean density of eigenfrequencies (modes) of the uncoupled subsystems, 
and is the energy stored in each subsystem. In a similar way, by recalling 
the average power ratio ( l28l) . it is worth noticing that the difference between 
the power entering two arbitrary elements of the chain can be written in the 
very general form 

p. T'C^.i) p. xihj) 

^ ~ ' ^^^^ 

valid for an arbitrary chain section of two elements. In f p2|) . Pi is the input 
power entering cavity i, and Pj is the input power entering cavity j. In equi- 
librium, the net power exchange between the to subsystems is established 
by the power dissipated in each cavity, and by the fraction of power leak- 
ing through the small (deterministic) channel. The scheme of the physical 
framework analyzed to derive (l32l) is reported in Fig. [71 By using ([5]) the 
loss factor \ and by assuming 

Pi P 

where L is a dissipation factor of the subsystem since ■* is the dissipated 
power within each cavity, it is found that 

' 27r \d. dJ ^ ' 
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cavity 



cavity 
di 



Figure 7: Arbitrary section of a chain: two cavities are excited by different 
input powers and they are weakly coupled. This simplified framework has 
been used to calculate the net power exchanged by the two cavities. A 
similar attack is adopted by the statistical energy analysis (SEA) approach 
in acoustics. 



where a coupling loss factor rjij 



= LT^'^^y27T is recognized, and 



(35) 



is the mean (mode) density of each subsystem. The assumption in fl33|) is 
related to the equipartition of energy. This results in f lM|) being formally 
similar to ( 131]) even in the case of autonomous (sub) systems coupled by an 
electrically short transmission line. Interestii igly, in deriving the t erms of 
(^^, a weak coupling argument has been u sed iGradoni et al.l (l2012dl ) that is 
identical to the one typically used in SEA Langleyl ( 1990l ). Methods based 
on SEA have been recently integrated with the f inite element r netho d (FEM) 
to tackle very complex vibro-acoustic scenarios iTanner et al. 



7. Conclusion and future perspective 

The general "dressed" impedance obtained in formulating the ROM for 
ports and terminals radiating inside complex electromagnetic cavities has 
been described. This statistical model accounts for short orbit deviations 
from the full chaotic mixing assumption. Numerical and experimental results 
validate the "dressing" of a universal fluctuating impedance matrix (depen- 
dent on a single loss parameter) with radiation impedance/admittance ma- 
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trices (modeling radiation of ports/apertures). A generalization of the RCM 
to the scenario of interconnected cavities for weak coupling conditions is 
presented, and the distribution of the power flowing through such a chain 
discussed for an increasing number of cavities. The analogy with Statistical 
Energy Analysis used in acoustics is pointed out by formulating a new model 
for the power exchanged by two wave chaotic cavities excited by different 
input powers. 
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